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ABSTRACT 

We present a model for the formation of massive (M> IOMq) stars through accretion- 
induced collisions in the cores of embedded dense stellar clusters. This model circum- 
vents the problem of accreting onto a star whose luminosity is sufficient to reverse the 
infall of gas. Instead, the central core of the cluster accretes from the surrounding 
gas, thereby decreasing its radius until collisions between individual components be- 
come significant. These components are, in general, intermediate-mass stars that have 
formed through accretion onto low-mass protostars. Once a sufficiently massive star 
has formed to expel the remaining gas, the cluster expands in accordance with this 
loss of mass, halting further collisions. This process implies a critical stellar density 
for the formation of massive stars, and a high rate of binaries formed by tidal capture. 
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1 INTRODUCTION 

Star formation on galactic and extragalactic scales is largely 
concerned with the formation of massive stars, due to their 
radiative, kinetic, and chemical feedback into the interstel- 
lar medium. However, the formation of massive stars is 
an, as yet, unsolved problem in astrophysics. Although 
low and intermediate mass stars (M < IOMq) are readily 
explained through gravitational collapse and subsequent ac- 
cretion (Palla & Stahler 1993), this model fails for high-mass 
stars. Radiation pressure, from a ~ IOMq stellar core, on 
the dust in the infalling gas, halts accretion and thus limits 
the mass (Yorke & Krugel 1977; Wolfire & Cassinelli 1987; 
Beech & Mitalas 1994). The aim of this paper is to explore 
an alternate model that circumvents this problem. 

A valuable clue to the formation mechanism of massive 
stars is found in their environments. Massive stars are rarely 
found in isolation. Rather, they are predominantly found 
in the central regions of rich stellar clusters (Zinnecker et 
al. 1993; Hillenbrand 1997). Furthermore, even the few run- 
away OB stars are best explained as having originated in 
the centres of dense systems (Clarke & Pringle 1992). 

Accretion in stellar clusters is likely to be an impor- 
tant mechanism in shaping the initial mass function (Zin- 
necker 1982; Larson 1992). Recent work on the origin of 
the stellar mass distribution has shown that competitive ac- 
cretion in clusters naturally results in the formation of the 
highest mass stars in the deepest parts of the cluster poten- 
tial, ie, the centre (Bonnell et al. 1997). The alternative, 
dynamical mass segregation, is unable to account for the lo- 
cation of the massive stars in such young clusters as the ONC 
(Bonnell & Davies 1997; Hillenbrand & Hartmann 1997). 



Although this model of competitive accretion cannot apply 
to very massive stars, where radiation pressure impedes the 
growth of stellar mass beyond 10 M@, it can qualitatively 
account for the observed trend of the more massive stars to 
be more centrally located. In the present paper, we combine 
the model of accretion onto low-mass stars in a protoclus- 
ter with the effects that such accretion has on the stellar 
dynamics of the cluster. We show that in this case stellar/ 
protostellar collisions become important in the densest part 
of the cluster which can lead to the build-up of very massive 
stars. 

In Section 2, we outline why forming massive stars is 
problematic. In Section 3, we discuss how the core of a clus- 
ter evolves during accretion. Section 4 investigates when col- 
lisions become important and Section 5 presents the models 
for the formation of massive stars. The question of binaries 
among the massive stars is addressed in Section 6. Section 
7 explores the implications of this process while our conclu- 
sions are presented in Section 8. 

2 THE PROBLEM 

There are a number of reasons why forming massive stars 
is problematic within our current understanding of star for- 
mation. Firstly, in isolated formation scenarios, there is 
the problem due to the difficulties of accreting onto mas- 
sive protostars. The spherical symmetric collapse of massive 
molecular clouds can be halted once the protostellar core has 
attained masses M > IOMq . This occurs due to the large lu- 
minosities of these stars, such that the radiation pressure on 
dust, which transfers its momentum to the gas, can halt the 
collapse and reverse the infall (Yorke & Krugel 1977; Wolfire 
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k. Cassinelli 1987; Yorke 1993). This limits the amount of 
mass that can be directly accumulated onto a massive star. 
For accretion to continue onto a massive protostellar core, 
the dust properties (abundances and sizes) have to be signif- 
icantly reduced from those found in the interstellar medium 
(Wolfire & Cassinelli 1987). 

This limit applies to spherically symmetric collapse and 
accretion. Generally there is some angular momentum in 
the parent cloud and this results in the formation of a pro- 
tostellar disc (eg, Yorke, Bodenheimer & Laughlin 1995). 
Accretion from this disc onto the massive protostar may be 
able to proceed if the disc is sufficiently thin, if not the disc 
will be destroyed (Yorke & Welz 1996). Unfortunately, it 
is presently unclear what disc thickness is necessary or ap- 
propriate. Still, to be able to form a massive star this way 
requires that all the gas collapses to form a disc before the 
protostar has gained a large fraction of its eventual mass. 
The problem with this is that collapse is non-homologous, 
with the central, low angular momentum regions collapsing 
faster. Plus, gravitational torques will ensure that the disc 
mass is always less than that of the central object, as they 
act on a timescale of order the disc's dynamical time, much 
shorter than the cloud's free-fall time, % (Bonnell 1994; 
Laughlin & Bodenheimer 1994). Thus, a central massive 
protostar should form while the majority of the cloud is still 
infalling, and this limits how much mass gathers in the disc 
and hence the efficiency of the process. 

Even if direct collapse could form massive stars, it is 
difficult to see how this process can account for the massive 
stars in the centre of the dense clusters such as the Trapez- 
ium in the Orion Nebula Cluster [ONC]. This is because 
the Jeans mass in the Trapezium protocluster cloud must 
have been small, perhaps as small as 0.3 M©(Zinnecker et 
al. 1993). To see this, we note the Jeans length (fragmenta- 
tion scale), 




must have been smaller then the average half separation of 
the stars that are now present in the cluster center ("touch- 
ing protostars" ) , where T is the temperature, p is the den- 
sity, R g is the gas constant, G is the gravitational constant, 
and fi is the mean molecular weight. From the current stel- 
lar number density (sa 2 — 4 x 10 4 stars pc~ 3 [McCaughrean 
& Stauffer 1994; Hillenbrand 1997]), we estimate i?j < 3000 
au. 

We can also estimate the mass density (« IOOMq now 
in stars) inside the current core radius (7? C orc ~ 0.05 pc, 
McCaughrean & Stauffer 1994) of the Trapezium, yielding 
~ 2 x 10 5 Mq pc -3 , equivalent to an initial smeared out gas 
density of p w 10~ 17 g cm~ 3 (or more, if the star formation 
efficiency, SFE, was less than 100 % and some gas was lost). 
With Rj and p as above, it follows then that the Jeans mass 
is small, Mj < 3M Q (assuming 100% SFE), where 




This implies a temperature of 90 K. A more realistic temper- 
ature of 20 K yields Mj < O.3M0. To raise the initial Jeans 
mass to 50 Mq, the mass of Ori 6 1 C, assuming the mass den- 
sity of the Trapezium cluster core, requires a Jeans length 



of 8000 au or 0.04 pc (larger than the present day star-star 
separations) and a temperature of T > 500 K. These values 
are clearly implausibly high. 

Thus, the stars that form in the centre of dense clus- 
ters should have initial masses much smaller than the me- 
dian mass in the final cluster. Their initial masses are un- 
likely to be larger than the minimum mass in the cluster. 
In the ONC, this corresponds to an initial Jeans mass of 
Mj ~ O.1M0. Furthermore, we know that the massive stars 
in the centres of young clusters had to form in situ as the 
clusters are generally not old enough to allow for dynamical 
mass segregation (Bonnell & Davies 1998). To obtain the 
observed high masses, subsequent processes including accre- 
tion onto the protostars, and, once M > 10M Q , collisions 
between the intermediate mass protostars, are required. 

The problem of forming a massive star through direct 
collapse and accretion also applies to scenarios where sig- 
nificant collisional buildup of gaseous clumps occurs before 
the collapse (eg Bastien 1981; Murray & Lin 1996). Instead, 
collisions need to occur after the objects have collapsed and 
attained stellar densities. Random encounters in star form- 
ing regions may lead to some protostellar collisions (eg Price 
& Podsiadlowski 1995) but will in general be very rare. In 
order for collisions to be significant, a high density environ- 
ment is required. 

3 ACCRETION AND CLUSTER DYNAMICS 

The presence of a large amount of gas in a young stellar 
cluster (e.g. Lada, Strom & Myers 1993) significantly af- 
fects the cluster's evolution. This gas, typically comprising 
1.5 to 10 times the mass in stars when the cluster is first 
detectable in the IR, is a remnant of the star formation pro- 
cess, and must have been significantly higher earlier on (ie 
100 % gas initially). Accretion of this gas by individual stars 
increases their masses, and, results in a spectrum of stellar 
masses even from an initially uniform distribution (Bonnell 
et al. 1997). 

Of direct significance here is the effects accretion has on 
the cluster dynamics. In the absence of large scale motions 
such as rotation, the gas velocities will be uncorrelated with 
those of the stars (especially in the cluster core). Thus, the 
accreted gas will, in general, have zero net momentum com- 
pared to the stars. Now, as each star accretes from the zero 
momentum gas, it conserves momentum and thus deceler- 
ates. This decreases the velocity dispersion in the cluster 
and hence decreases its kinetic energy. Additionally, the 
increase in the stellar mass decreases the cluster's potential 
energy. Thus, the combined effect of accretion in a cluster is 
to decrease the total energy, and thus make it more bound. 
This causes the cluster to shrink to compensate for its new, 
lower, total energy. 

For the purposes of this paper, we consider the effects 
of accretion on the dynamics of the core of a young stel- 
lar cluster. The core is modelled as a collection of N stars 
within a core radius Rcore, total stellar mass M s tars, and 
total core mass M CO ro, M CO re = M sta rs + M gas , where M gas 
is the mass in gas. The individual components of the core 
are not modelled. Instead, the core's evolution is followed 
by considering its global properties and how they evolve un- 
der accretion. We consider the core as the central region 
with a shallow radial density distribution, embedded within 
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a larger, steeper (such as an R~ 2 ) distribution. Thus, we 
assume we can neglect the contribution to the potential and 
kinetic energies from the surrounding cluster. We also as- 
sume that the number of stars in the core is fixed (no stars 
enter or leave) and that the gas can only enter (and not 
leave) the core until a massive star is formed. 

In the absence of accretion, the core will be in virial 
equilibrium such that the core radius is 

GA/corc 



R, 



(3) 



disp 



where v disp is the velocity dispersion in the core. The pri- 
mary effect of accretion is to perturb the core from this 
equilibrium by decreasing the total energy. The core then 
adapts by shrinking, until the velocity dispersion increases 
sufficiently to reestablish virial equilibrium with a new R,corc 
and Vdisp- These quantaties can be calculated by first con- 
sidering how the total energy decreases as the mass is added 
and secondly by determining what core radius Rcore (and 
vjisp) are appropriate for a virialised system at the new, 
lower, total energy. 

In order to calculate how the cluster core shrinks with 
accretion, we consider the total energy of the stars. The ki- 
netic energy of the stellar component of total mass M s t ars j is 
-Ekin = l/2M s tars fdispi anc ^ ^he potential energy of the stars 
from both the stellar and gaseous components is E po t = 
-GM sta rs (A/stars + M gas )/-Rcoro. To see the results of ac- 
cretion more directly, we express the stars' kinetic energy in 
terms of a characteristic stellar momentum, p, defined from 
the velocity dispersion as, p — M s tars«disp- The total energy 
of the stars, in a cluster core of total mass M CO re, is then 
given by: 

„2 



P 



GM COIC M s tnrs 



(4) 



2-Mstars ^corc 

From equation 4, we see that adding mass to the core while 
conserving momentum decreases the total energy. The stars 
become more bound (— GM core Mstars/-Rcore) as they have 
a larger mass. Similarly, their kinetic energy (and velocity 
dispersion) decreases due to the conservation of momentum, 
p, (v disp oc M~l IS , E kiD = p 2 /2M st ars oc M^J. This 
new lower total energy, and lower kinetic energy means that 
the core is now out of virial equilibrium. It then shrinks 
while conserving energy (increasing the kinetic energy and 
the Vdisp) until virial equilibium is restored. Once the cluster 
core revirialises with total energy _E v ; r (the decreased -Btot 
after accretion), its radius is given by 

— GAfcore Af s tars 



-Rcorc- 

^J_/y lr 

and the new velocity dispersion is 



^disp 



GM C . 



Rc 



(5) 



(6) 



To follow the evolution of the cluster we use equation 4 
to calculate the destabilisation of the cluster core due to 
the accretion, and then equation 5 to calculate the new core 
radius once it has revirialised. The two remaining unknowns 
that need to be considered are the timescale on which this 
happens and the mass of the gas contained in the cluster 
core. The minimum timescale for any process considering 
the dynamics of a group of stars is the crossing time, icross = 



J?/ Vdisp, where R is a characteristic radius of the stellar 
distribution. 

The mass of gas (M gas = M CO rc — Af s tars) contained 
within the core is uncertain and depends on the cluster 
dynamics and how the gas reacts to the deepening poten- 
tial. There are three main possibilities. The first is that 
the accreted gas comes primarily from outside the cluster 
core, such that the gas mass-fraction in the core is constant. 
The second is that the accreted gas is completely contained 
within the core and thus the core mass is constant while the 
gas mass decreases during accretion. The last possibility is 
that the gas infall onto the core is greater than the accretion 
onto the stars and thus the fraction of the core mass in gas 
increases. Although it is unclear which of these scenarios is 
correct, we can make a qualitative estimate based the fact 
that the gas infall time should be comparable to the dynam- 
ical, or crossing, time (icross = R/ Vdisp = y/Rf/GM). In 
the near-uniform density core of the cluster, M oc R 3 , the 
crossing time does not change drastically with radius so the 
gas has time to catch up so the gas can fall in at the same 
rate as the core's radius decreases. In the surrounding areas 
with a steeper (R~ 2 ) density profile, M oc R, the crossing 
time increases rapidly with increasing R, and gas that gets 
left behind cannot catch up to the core. In reality, the core 
can only shrink with accretion so it will wait for the gas 
from outlying areas to fall into it, but it does mean that the 
core is never overwhelmed by gas from outside, it will adapt 
as quickly as the gas is fed in. We assume that no gas can 
leave the core; that no dominant outflows capable of remov- 
ing a significant fraction of the gas occur before the massive 
stars are formed. Thus, any gas which is initially in the core 
will remain in the core until accreted. Furthermore, as the 
accretion-induced contraction of the cluster requires several 
initial crossing times (see below) , significant amounts of gas 
will fall in from outside the core. Therefore, for the rest 
of this paper, we generally assume that the gas primarily 
comes from outside the core, and that the core gas mass is 
proportional to the stellar mass. We do, however, comment 
on how the different scenarios affect the results. 

The rate at which the cluster core shrinks with added 
mass depends on how the accretion time, t a cc, compares to 
the crossing time, t cr oss, as the core will revirialise on the 
order of its crossing time. For example, if the mass is added 
on timescales t acc << tcross, no matter how much mass is 
added, the cluster core will revirialise at no less than half 
its original radius. 



3.1 Slow Accretion 

For slow mass accretion, t acc >> tcross, the core adapts 
practically instantaneously. In this case, the change in ra- 
dius as a function of accreted mass, can be derived by dif- 
ferentiating both equations 4 and 5 with respect to mass, 
remembering that the change in -Btot with mass occurs at a 
constant radius whereas the change in E v - U reflects changes 
in both mass and radius, and that momentum is conserved 
during accretion. When the core revirialises, the change in 



i? v ir and Etot are equal: 



dE t , 



dM si 



dE v 



dM s - 



(7) 
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Using the fact that p /M st ars = GMcoreMstarsZ-Rcorc when 

virialised, and assuming M CO rc oc M s tars, the core radius can 
be shown to decrease due to accretion as: 



-Rcore OC M , 



(8) 



Alternatively, if the core mass is assumed to be constant 
and the accreted gas is initially fully contained within the 
cluster core, then M CO re = -Mstars + -Mgas is constant and 
the cluster shrinks due to the loss of specific kinetic energy. 
In this case 



Rc 



oc M~ 



(9) 



In the third scenario where the infall of gas is such that 
it more than compensates for the accretion, then M CO rc oc 
Mgas, and the cluster core shrinks primarily due to the de- 
crease in potential energy. In this case 



7?core OC M~ 



(10) 



The decrease in Rcoro in this case will then be from a com- 
bination of the gas infall rate (equation 10) and the stellar 
mass accretion rate (equation 8). 

In all cases, it is clear that even small amounts of ac- 
creted mass can significantly decrease the core radius, due to 
the decrease in total energy (equation 4) , and this increases 
the stellar density. This is important because it dramatically 
increases the rate of collisions (see §4 below). 

3.2 Fast Accretion 

For moderate and high mass accretion rates, face Across, the 
cluster core does not revirialise immediately, and significant 
amounts of mass are added before the core can adjust. In 
order to quantify how this changes the evolution of i? C orc, 
we performed several direct nbody simulations (using a tree- 
code derived from the SPH/N-body code of Bate, Bonnell 
& Price 1995) of clusters of 100 stars within a radius of 0.1 
pc undergoing constant accretion. These simulations were 
compared with a toy model of the cluster core's evolution un- 
der the influence of accretion. This model calculates, first, 
how the core shrinks as the total energy decreases under 
accretion, assuming the mass is added at a constant core 
radius (equation 4). The amount of mass added is equal 
to the accretion rate multiplied by the timescale for the 
process (the revirialisation timescale). Secondly, the clus- 
ter then revirialises at a radius corresponding to the new, 
lower, total energy (equation 5). The timescale for the re- 
virialisation is a free parameter, on the order of a crossing 
time. By comparing the N-body simulation with this toy 
model, we determined that the cluster can effectively re- 
virialise every 0.75 tcross (Figure 1). In both cases, the ac- 
creted gas is assumed to come from outside the core such 
that Mcore = M s tars (no unaccreted gas in core). The two 
cases are well matched after w ltcross, showing that the sim- 
ple model for the cluster's evolution under mass accretion 
is a good representation of the dynamics involved. Note 
that once the cluster has shrunk, the crossing time is much 
less than the accretion timescale and the variation in cluster 
radius follows equation 8. It is also worth noting that al- 
though the velocity dispersion initially decreases under ac- 
cretion (due to the conservation of momentum), the sub- 
sequent revirialisation will increase the velocity dispersion 




Mass 

Figure 1. The half-mass radius, in units of the initial half-mass 
radius, versus total stellar mass, in units of the initial cluster 
mass, for the nbody simulations (solid lines) and the toy cluster 
model (dotted lines). The lines are for (from left to right) mean 
stellar accretion rates of 1 X 10~ 7 , 5 X 10~ 7 , 2 X 10~ 6 , 5 X 10~ 6 
and 1 X 10~ 5 MQvr" 1 . 

as potential energy is transferred into kinetic energy (equa- 
tion 6). For isolated clusters, the comparisons are only valid 
for approximately a relaxation time (several crossing times 
for 100 stars), as an isolated cluster of 100 stars starts to 
dissolve (due to the formation of a hard binary which can 
absorb the total energy of the cluster). However, if the 100 
stars are embedded in a much larger system, it does not dis- 
solve, but simply exchanges members with the surrounding 
cluster. 

Thus, for most accretion rates, the rate at which the 
cluster (and core) shrinks depends not only on the accretion 
itself but also on the crossing time. This is important as 
it decreases the differences between the rates of shrinkage 
(equation 8, 9 and 10) for the various gas dynamics sce- 
narios. The crossing time depends on the amount of mass 
in the cluster core, such that the crossing time is signif- 
icantly shorter when all the accreted gas is initially con- 
tained within the core as the velocity dispersion is larger. 
Thus, even though the core does not shrink as much with 
accreted mass in this case (equation 9), the fact that it can 
adapt much more quickly balances this difference and the 
core actually evolves in a similar manner. In the following 
sections we assume that the accreted gas primarily comes 
from outside the cluster such that M CO ro oc M s tars- 



4 TIMESCALES: 
ACCRETION 



COLLISIONS VERSUS 



In order to ascertain how and when collisions become im- 
portant, we need to quantify and compare the timescales 
for collisions and accretion. The accretion timescale, t a cc, 
defined as the time it takes to double a star's mass is simply 
given by 

tacc = ^, (11) 

M* 
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where M* is the star's mass and M* is the stellar mass accre- 
tion rate. The collisional timescale has a more complicated 
form (see Binney & Tremaine 1987 for a derivation), de- 
pending on the stellar velocity dispersion, v^isp, the stellar 
density n, and the radius at which mergers occur, i?* : 



2xl0 5 4xl0 5 



1 

^coll 



= 16y^nv disp Rl ^1 + 



(12) 



The last term incorporates the effects of gravitational focus- 
ing in the cluster. Gravitational focussing has an important 
effect on the collision rate for the cases considered here. 

For collisions to be at all relevant, the cluster core must 
be sufficiently dense. Accretion onto the core of a cluster 
decreases its radius and increases its density until collisions 
eventually become significant. Figure 2 illustrates the evolu- 
tion of the accretion and collisional timescales for a cluster 
of 100 stars initially contained within a radius of 0.1 pc. 
The cluster members are all initially 0.1 M© and accrete at 
a rate of 2 x 10~ 6 M©yr _1 . The collisional timescale is cal- 
culated assuming stellar merger radii of 0.1 au and 0.01 au 
(R* « 21.5 and 2.15 R©). We assume that M CO ro oc M s tars 
and neglect the mass of any gas in the core. The main effects 
of a non-zero gas mass are to decrease the rate of shrinkage 
with accretion and to decrease the crossing time. The net 
effect is little different from Figure 2. 

Figure 2 illustrates the relevant timescales for doubling 
a star's mass (tmass)- The accretion timescale remains ap- 
proximately constant at a few x 10 5 years, while the collision 
timescale decreases from > 10 8 years to < 10 5 years. Note 
that the initial cluster crossing time is tcross = 1.5 x 10 
years and that in both cases, the t co u becomes smaller than 
tacc in < 3t C ross, or < 4.5 x 10 5 years. Thus, collisions be- 
come the dominant mode of increasing a star's mass in the 
cluster core within 5 x 10 5 years. 

The question of what the stellar merger radius should 
be is not easy to answer as the protostars considered are 
still in their pre-main-sequence phase, accreting from the 
surrounding material and their accretion discs. They are 
in general significantly larger and fluffier than their main 
sequence counterparts. Pre-main sequence stars, of various 
masses, less than 10 6 years old, most probably have radii 
in the range 5 to 20 R (Palla & Stahler 1993; Beech & 
Mitalas 1994). Furthermore, the presence of an accreting 
envelope and disc aids in the capture and merger process 
(eg Clarke & Pringle 1993; Hall, Clarke & Pringle 1996), 
thus giving the stars an effectively larger radius. For the 
remainder of this paper, we assume a merger radius of 0.1 
au (21.5 R©). ft is worth noting that a smaller merger 
radius still allows for the collisional buildup of massive stars. 
However, it does imply a larger mean stellar mass at the 
time when collisions become important as it allows for more 
accretion. 



5 COLLISIONS AND MASSIVE STARS 



In the previous sections, we have seen how the central cores 
of stellar clusters undergoing accretion eventually shrink 
sufficiently that collisions dominate the increase in stellar 
masses. In this section we use the same toy model described 
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Figure 2. The timescale to double a star's mass either by colli- 
sions (solid and long-dashed lines) or through accretion (dotted 
line) is plotted versus time in units of the initial crossing time, 
tcross = 1.5 X 10 5 years, (and in years) for a cluster of 100 stars 
with an initial radius of 0.1 pc. The accretion timescale assumes 
an accretion rate of 2 X 10 — ® M0yr — The merger radius for 
collisions is taken to be 0.1 (solid line) and 0.01 au (long-dashed 
line). The cluster radius in units of the initial radius is also plot- 
ted (short-dashed line). 

above but modified to include collisions and competitive ac- 
cretion (Zinnecker 1982). Accretion is allowed onto all stars 
in proportion to their contribution to the gravitational po- 
tential, ie 



Mi, 



Ml 



r^dust- 



(13) 



This is equivalent to each star accreting with a Bondi-Hoyle 
type accretion rate (eg Ruffert 1996) assuming that the ac- 
creting material is at a uniform density and each star is 
moving through the medium with the same velocity. The 
collisions occur on the collisional timescale of each star, t co n 
(equation 12), and the stars are assumed to merge when 
they collide. The collisions are assumed to be 100 % effi- 
cient in that no mass is lost. In reality, some small amounts 
of mass will be lost during collisions (eg Davies et al. 1993). 
Mass loss rates of up to 25 % of the stellar masses were 
attempted with no significant differences in the results. Ac- 
cretion is not allowed onto stars more massive than 10 M© 
so that stars more massive than this need to form through 
collisions. The cluster evolution is followed, allowing for re- 
virialisation every 0.75t C ross, and recalculating t co u for each 
star. 

The number of stars is only affected by collisions. No 
stars are allowed to leave or enter the cluster core. In real- 
ity some interchange between the core and the surrounding 
cluster is expected, introducing some lower-mass stars and 
removing some of the high-mass products. Ejection of stars 
from the core require an interaction with a hard binary in 
order to give a kick velocity greater than the velocity dis- 
persion. The finite and, including circumstellar material, 
relatively large size of the young stars limits such ejections 
as the interaction will often lead to a merged object (Davies, 
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Figure 3. The four panels represent the toy-model calculations of the central core of a cluster containing 100 stars initially within 0.1 
pc, evolving under mean stellar accretion rates from 10 — ''tolO - 4 M0yr — The stellar merger radius is taken to be 0.1 au. The top 
left panel shows the time, in years, at the end (dashed line) of the evolution (either 10 6 years or when only 2 stars are left), and the 
time when a 50 Mq star is first formed (heavy solid line). The top-right panel shows the maximum stellar mass (in Mq) for the same 
cases plus the median stellar mass at the end of the evolution (light solid line). The bottom left panel shows the number of collisions 
that occur for the two cases while the radius (in pc) of the cluster core is shown in the bottom right panel. 



Benz & Hills 1994). Furthermore, the presence of the gas 
in the cluster should decrease such ejections and, in gen- 
eral, it would be preferentially the less massive stars that 
are lost from the core. If many interchanges occur, then the 
introduction of low-mass stars will increase the number of 
collisions required to form a massive star. In the extreme 
case, where many more stars leave the core then enter it, an 
initially higher number of stars in the core may be necessary. 
These additional complications are beyond the scope of the 
present study but should not overtly change the results. 

The models are followed for 10 6 years or until only two 
stars remain (a rare occurrence). Models with varying num- 
bers of stars in the cluster core and varying core radii were 
investigated. In general, cluster cores with t cr oss signifi- 
cantly less than 10 6 years (the approximate age at which a 
complete mass spectrum is observed in the ONC) are the 
ones in which significant collisional buildup of massive stars 
occurs. 



An illustrative case of a cluster core with 100 stars in 
0.1 pc, initial stellar masses of OTM© and a stellar merger 
radius of 0.1 au, evolving under various mean stellar accre- 



tion rates is shown in Figure 3. This figure shows the effects 
of accretion-induced collisions in forming massive stars. The 
top-left panel shows the period of time required to form a 
star of 50 Mq, and the time at which only 2 remain of the 
original 100 protostars (or 10 6 years). It is worth noting that 
the evolution requires several x 10 5 to 10 6 years to form a 50 
Mq star. Thus, massive stars in these clusters are expected 
to be considerably younger than the mean stellar age, and 
in some cases, could be up to 10 6 years younger than the 
oldest low-mass stars. 

The bottom-left panel of Fig. 3 shows the number of 
collisions that occur and the bottom-right panel shows the 
cluster's core radius for the same cases. For small mean stel- 
lar accretion rates, few collisions occur as the cluster does 
not shrink sufficiently within 10 6 years. For large accretion 
rates, the mass is added very quickly, face << Across, and 
by the time the cluster revirialises, the mean stellar mass 
is greater than 10 Mq, halting the accretion. In this case 
the cluster core does not shrink very much (half the initial 
core radius) and hence collisions are rare. Maximal shrink- 
age occurs when i a cc ~ Across ~ 10 6 years. Once collisions 
dominate, the cluster evolves due to the collisions and has 
no time to shrink further. This occurs at slightly larger radii 
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Figure 4. Same as Figure 3 except for a cluster core of 100 stars initially within 0.15 pc and with a gas fraction of 50 % contributing to 
the gravitational potential. The gas is removed once a star with M >50Mq has formed. The removal occurs over several crossing times. 



for larger accretion rates and the mean stellar mass is then 
larger. 

The top-right panel of Fig. 3 shows the maximum mass 
obtained at either 10 6 years, or when only 2 stars are left, 
or when a collision has formed a star of mass > 50 Mq. It 
is clear from this panel that the maximum mass depends on 
both the accretion rate and the number of collisions. The 
maximum mass is greater than 50 Mq for accretion rates 
between 5 x 10~ 7 and 5 x 10~ 5 Mq yr^ 1 . The median mass 
at the time that a 50 Mq star is first formed is also shown 
in Fig. 3. The median mass typically lies between a few and 
10 Mq. It is worth noting that for accretion rates > 10~ 5 
Mq yr _1 , accretion alone will produce a 10 Mq star in 10 6 
years. 

As noted above, the exact value of the stellar merger ra- 
dius does not have a large affect on the results. The above 
calculations were also performed with a stellar merger radius 
of 5 Rq and gave qualitatively similar results, the main dif- 
ferences being the amount the cluster shrinks as the smaller 
stellar merger radius necessitates a higher stellar density for 
collisions, and the mean stellar mass which is larger in the 
smaller stellar merger radius case. 

The maximum mass in Fig. 3 sometimes reaches well 
in excess of 100 Mq. This is clearly not physical and is 
due to the artificial termination points for the evolution (2 
stars left, or 10 6 years). A more realistic termination point 
requires some process that stops the collisions. Up to now, 



we have ignored the contribution to the potential of the 
gas being accreted. The gas component in young clusters 
typically comprises 50 to 90 % of the total mass in the cluster 
(Lada et al. 1993). The exception to this is the Trapezium 
where the most massive star, 9 C, has presumably removed 
all the gas through its ionising photons and wind. This offers 
an appealing process of halting further collisions. 

In order to include the effect of the gas on the potential, 
we performed several calculations with either a gas contri- 
bution equal to that of the stars or a constant core mass and 
thus an ever decreasing gas mass. In both cases, the gas is 
presumed to be removed over several crossing times once a 
star of sufficient mass (50 Mq in this example) is formed. 
The gas is expelled due to the ionising photons and wind 
from the massive star. 

Firstly, we consider the effect of the gas when it con- 
tributes equally with the stellar mass. Figure 4 shows the 
evolution in this case for a cluster core of 100 stars initially 
within a radius of 0.15 pc (the increased radius balances the 
added mass component of the gas, so that the crossing time 
is similar). The evolution is very similar to the previous 
one shown in Figure 3 until a 50 Mgstar is formed. At this 
point, the gas is removed and the cluster relaxes, increasing 
its radius, to account for the loss of potential energy. The 
larger radius, and hence smaller stellar density, effectively 
eliminates any further collisions. This impedes the forma- 
tion of the very massive stars (M > IOOMq). 
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Figure 5. Same as Figure 3 except for a cluster core of 100 stars initially within 0.15 pc and with a constant core mass. The gas is 
initially contained within the core and has mass equal to that required for the accretion to be sustained for a period of 2 X 10 6 years. 
The gas is removed once a star with M > 50Mq has formed. The removal occurs over several crossing times. 



Alternatively, if the gas is completely contained within 
the core initially, and the mass of the core is constant, then 
the core radius does not shrink as rapidly with accreted mass 
(Rcore oc M~ t ^ rs ), but this is somewhat offset by the much 
shorter crossing time as the velocity dispersion is higher due 
to the larger core mass. Such an evolution is illustrated in 
Figure 5. The initial gas mass is set to twice that required 
for the desired accretion rate over 10 6 years. Thus, the star 
formation efficiency increases from near zero to a maximum 
of 50 %. The minimum cluster radius is well represented by 
equation 9: -Rcorc/-Ri n it = {M* x t/Af*; n j t )~ 2 . The main 
difference between the cluster evolution in this case and in 
the previous ones is that, although the cluster only shrinks 
as , due to the large amounts of gas included in the 

cluster core, the initial crossing time is much shorter. Thus, 
it can adapt very quickly and there is no maximum value of 
the mean stellar accretion rate for which significant shrink- 
age occurs. Any large rate (M*>10~ 6 MQyr^ 1 ) is sufficient 
to drive the collisional buildup of massive stars. The large 
amounts of gas included in the cluster core result in the core 
being unbound and expanding once the gas is removed by a 
newly formed massive star (M* > 50Mq). In this case, the 
massive stars quickly expand into the rest of the cluster. 



6 BINARY SYSTEMS 

So far, we have neglected binary systems and their possible 
role in the evolution of the cluster core and in the formation 
of massive stars. Binary systems are commonly employed 
as an energy source in globular clusters to stave off core- 
collapse. The binaries inject energy through binary-single 
interactions transforming their potential energy into the ki- 
netic energy of the cluster. 

In the scenario considered here, the cluster is shrink- 
ing due to the accretion. The rate of shrinkage depends on 
the relation between the accretion timescale and the cross- 
ing time. In the extreme case, t acc » icross, the shrinkage 
goes as i?core oc MZgrg- This is the same rate as is found 
for a binary system that is accreting zero-angular momen- 
tum matter (Smith, Bonnell & Bate 1997). Thus, unless a 
primordial binary is initially capable of affecting the cluster 
dynamics, the subsequent accretion and shrinkage should 
not make it an important energy source for the cluster core. 

We do know, however, that many massive stars are in 
close binary systems (Garmany, Conti & Massey 1980, Ma- 
son et al. 1996) so that any theory for their formation should 
be able to explain why they are in such systems. Their mul- 
tiplicity can easily be explained in the context of the model 
presented here where massive star formation is due to col- 
lisions. For collisions to be important, many more interac- 
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tions will occur at slightly greater periastron separations. 
These interactions can result in the formation of binary sys- 
tems through tidal capture (eg Mardling 1996). Tidal cap- 
ture can produce binaries at a few stellar radii but will not 
form very wide systems. Furthermore, the next interaction 
can result in a loss of kinetic energy from the binary and 
thus help the stellar merger rate. In any case, the end re- 
sult is likely to be of a binary system comprising two of 
the remaining (proto)stars. The resultant binary is likely to 
comprise two relatively massive stars, as they both should 
have accreted a significant amount of mass, and/or under- 
gone collisional mass buildup, before capture. The number 
of binary systems should be large amongst the massive stars 
because the tidal capture rate must exceed that for colli- 
sions. 

Lastly, binary systems among the massive stars can pro- 
vide a test of this mechanism for the formation of massive 
stars. The maximum separation of such systems should not, 
in general, be larger than the hard/soft boundary typical for 
the stellar density when the collisions occurred. Except for 
the possibility of binary formation at the very end of the col- 
lisional buildup of massive stars, all binaries of separations 
greater than this critical, density-dependent, separation are 
broken up by encounters within the core. 

7 IMPLICATIONS 

There are several implications that can be gleaned from 
the models presented here. The accretion-induced collisions 
mechanism requires that the massive stars form in the centre 
of a large stellar cluster. Thus massive stars should never 
form in isolation, and those that are found in isolation need 
to have been ejected from such dense systems. 

Accretion-induced collisions can, under reasonable ac- 
cretion rates, form massive stars in less than 10 6 years. The 
process does require relatively dense initial conditions, in 
order to be able to accrete and shrink sufficiently quickly 
so that collisions are frequent. This is equivalent to requir- 
ing tcross ~ iacc << i C nd where t cnd « 10 6 years or some 
similar maximum time. For initial masses « O.IMq, this 
implies minimum stellar densities of >10 4 stars pc -3 . Such 
densities are typical only of the central regions of large clus- 
ters such as the Orion Nebula cluster, where these massive 
stars are indeed found (McCaughrean & Stauffer 1994; Hil- 
lenbrand 1997). 

The time required for accretion-induced collisions to 
form a massive star may imply a considerable age difference 
between the older low-mass stars and the younger high-mass 
stars. This age difference is at least several xlO 5 years and 
can be as much as 10 6 years, depending on the initial cluster 
density and the accretion rate. Thus stars such as (^C in the 
ONC may have only recently "turned on" and started to af- 
fect the gas in the cluster and in discs around the low-mass 
stars. The Trapezium is then the remnant of the central 
dense core where the collisions occurred and which has re- 
cently expanded to its present configuration due to the gas 
expulsion. 

The collisional buildup of massive stars also has im- 
plications for the intermediate-mass stars. It is these stars 
that, having attained their mass through accretion in the 
cluster core, are required to collide to form the massive 
stars. Thus, the formation of the massive stars implies a 



removal of intermediate mass stars. Systems with massive 
stars could then demonstrate a relative lack of intermediate- 
mass stars. There is a tantalising hint of the relative paucity 
of intermediate-mass stars in the ONC (Hillenbrand 1997) 
and in a few other clusters (Wilner & Lada 1991) although 
it could be due simply to statistical fluctuations. 

Finally, since, in this model, massive (m > WMq) and 
intermediate mass (m < WMq) stars are formed by differ- 
ent mechanisms, the resultant mass distribution may display 
some structure at « 10M© (not a single power-law IMF). 

8 CONCLUSIONS 

Accretion-induced collisions in the core of a dense, young 
stellar cluster can form massive stars. This circumvents 
the problem of accreting directly onto stars with masses 
> WMq where radiation pressure can halt the infall. Gas 
accretion onto individual members of the core of a stellar 
cluster forces the core to shrink. If the accretion timescale 
is comparable to the initial crossing time, the core can shrink 
sufficiently that collisions become significant. The collisions 
involve intermediate mass stars that are formed by accre- 
tion onto initially lower mass protostars. The actual stellar 
merger radius does not overtly affect this process as the 
cluster core will shrink due to the accretion until collisions 
become dominant. The collisions are most probably halted 
by the ejection of the gas contained in the core once a suffi- 
ciently massive star is formed (« 50M©). 

This mechanism for forming massive stars implies that 
they need to form in the centre of rich, dense young stellar 
clusters. They could appear significantly younger than the 
mean age of the low-mass stars. There is also the potential 
that the collisions will deplete the number of intermediate 
mass stars as it is they that will collide to form the massive 
stars. Close binary systems should be common amongst the 
massive stars, formed through tidal capture of stars which 
have close interactions but do not collide. 
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